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ABSTRACT 

The atmospheres of hot Jupiters and other strongly-forced exoplanets are 
susceptible to a thermal instability in the presence of ohmic dissipation, weak 
magnetic drag and strong winds. The instability occurs in radiatively-dominated 
atmospheric regions when the ohmic dissipation rate increases with temperature 
faster than the radiative (cooling) rate. The instability domain covers a specific 
range of atmospheric pressures and temperatures, typically P ~ 3-300 mbar and 
T ~ 1500 — 2500 K for hot Jupiters, which makes it a candidate mechanism to 
explain the dayside thermal "inversions" inferred for a number of such exoplanets. 
The instability is suppressed by high levels of non-thermal photoionization, in 
possible agreement with a recently established observational trend. We highlight 
several shortcomings of the instability treatment presented here. Understanding 
the emergence and outcome of the instability, which should result in locally 
hotter atmospheres with stronger levels of drag, will require global non-linear 
atmospheric models with adequate MHD prescriptions. 



1. Introduction 

A variety of gaseous exoplanets with strongly-forced atmospheres have been discov- 
ered by astronomers (e.g., Charbonneau 2009). These exoplanets, exemplified by the hot 
Jupiter class, receive extreme levels of irradiation from their stellar host and likely expe- 
rience permanent day- and night-side forcing conditions, from being tidally-locked on very 
compact orbits. Observationally, two of the most interesting trends emerging from studies 
of well-characterized, transiting hot Jupiters have been the tendency for many such planets 
to exhibit radius inflation, well above expectations from standard planetary cooling models, 
and the inference that thermal excesses (inversions) are present on the dayside of some of 
the strongly-irradiated planets (see Deming & Seager 2009; Baraffe et al. 2010; Burrows & 
Orton 2010; Winn 2010 for reviews). 
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It has recently been proposed that hot, strongly-forced exoplanet atmospheres are the 
site of significant magnetic induction when fast, weakly- ionized atmospheric winds cross 
the pre-existing planetary magnetic field. This induction takes the form of magnetic drag 
acting to brake atmospheric winds and associated ohmic heating in the planetary atmosphere 
and interior (Batygin & Stevenson 2010; Perna, Menou & Rauscher 2010a,b; Menou 2012). 
While other mechanisms have been proposed, ohmic heating is currently one of the leading 
contenders to explain the inflated radii of hot Jup iters (Batygin et al. 2011; Laughlin et al. 
2011; Menou 2012; Wu & Lithwick 2012). 

In this letter, we explore the possibility that thermal inversions in hot Jupiter atmo- 
spheres are not caused by extra absorption of stellar light at altitude, as they have tradi- 
tionally been interpreted (Hubeny et al. 2003; Fortney et al. 2006, 2008; Burrows et al. 
2008; Madhusudhan & Seager 2010), but instead result from the same induction mechanism 
that may explain radius inflation. In this alternative interpretation, thermal inversions have 
their origin in a thermo-resistive instability that affects radiatively-dominated atmospheric 
regions under specific conditions of weak magnetic drag, strong ohmic dissipation and fast 
winds. 



Let us first establish the existence of a thermo-resistive instability under physical con- 
ditions relevant to hot exoplanet atmospheres. 

Energy Equation - The energy equation satisfied by an atmosphere experiencing 
ohmic dissipation can be written 



where p is the mass density, C p is the specific heat at constant pressure and Qohm is the 
volumic ohmic dissipation rate. In the last expression, the divergence of the net radiative 
flux, F rad , has been rewritten as a net volumic radiative heating rate. 

In radiatively-dominated regions of an atmosphere, which is typically the case at and 
above infrared photospheric levels in hot Jupiter atmospheres (Seager et al. 2005, Cowan & 
Agol 2011, Menou 2012; Perna et al. 2012), the advective term (v • VT) can be neglected 
by comparison to the net radiative term (Q ra d)- Furthermore, near radiative equilibrium, 
this net radiative term can expanded into a first-order Taylor series around the local radia- 
tive equilibrium temperature, T eq (Goody & Yung 1989), so that the energy equation in a 
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radiatively-dominated atmosphere reduces to 



"■L eq , ^johm /-|\ 

<9t T rad pC p ' 

where r ra( j is the time to relax to radiative equilibrium in this so-called Newtonian ap- 
proximation. Ohmic dissipation provides a strictly positive contribution, so that in steady- 
state (dT/dt = 0), one expects thermal balance through net radiative cooling, at an atmo- 
spheric temperature T > T eq . This defines the basic radiative-ohmic state on which we wish 
to perform a stability analysis. 

Magnetic Induction Formalism - We use the steady-state axisymmetric magnetic 
induction framework introduced by Liu et al. (2008) and further discussed in the hot Jupiter 
context by Batygin & Stevenson (2010), Perna et al. (2010a,b) and Menou (2012). In this 
framework, a weakly-ionized zonal wind with velocity scale flows across the planetary 
magnetic field of surface strength B. This induces an additional magnetic field in the atmo- 
sphere and, in steady-state, this magnetic induction is balanced by resistive diffusion and 
associated ohmic dissipation in the atmosphere. 

Lacking detailed models for how electric currents flow in the generally inhomogeneous 
atmospheres of hot Jupiters and other hot exoplanets, we evaluate the heating term associ- 
ated with ohmic dissipation in Eq. (JT|) simply as 



H ^ = Qohm = V <f> B = Y± (2) 
pC p 4:TCT]pC p 

where rj is the electric resistivity and 7d rag = 4irrip/B 2 is the corresponding mag- 
netic drag time. The local electric resistivity of the atmospheric gas is evaluated as rj = 
230\/T /x e cm 2 s _1 , where x e is the free electron ionization fraction (e.g., Menou 2012). This 
expression for Q hm, which emerges from simple dimensional analysis of the induction equa- 
tion (Menou 2012), is consistent with more detailed calculations which find that the bulk 
of the ohmic dissipation is located in the active weather layer where induction occurs (e.g., 
Batygin & Stevenson 2010; Perna et al. 2010b; Wu & Lithwick 2012). The exponential 
dependence of 77 and Qohm with temperature when thermal ionization dominates is at the 
origin of the thermo-resistive instability. 

Linear Instability Criterion - Let us denote by C ra( j the first, radiative relaxation 
term on the right-hand side of Eq. <^). Radiative-ohmic equilibrium is satisfied when, in 
Eq. ([I]), the heating term H ohra (defined in Eq. [2]) exactly balances the cooling term C ra( j 
(with T > T eq ). This radiative-ohmic thermal equilibrium may be unstable, however, if 
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dT dT r rad ' 



so that a slight temperature increment would lead to a greater increase in ohmic heat- 
ing than in net radiative cooling. In the limit of weak drag, V<p changes so weakly with 
temperature, despite the varying magnetic drag (e.g., Menou 2012 or Eq. [3] below), that 

dH ohm ^ dH ohm drj _ H ohm drj 
dT drj dT rj dT 

To obtain a result in analytic form, we use here a simplified expression for the ionization 
fraction, which only accounts for potassium ionization (Balbus & Hawley 2000), 

Xp ^ = 6.47xW^V /2 f TV/4 
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'2.4 x 10 15 \ 1/2 exp(-25188/T) 
"nl J 1.15 x 10 



where n e and n n are the number densities of electrons and of neutrals, respectively (in cm -3 ), 
ax = 10~ 7 is the potassium abundance (assumed to be solar), and T is the temperature in 
K. With this expression for x e , we obtain 

= # ohm (2.5 x lO-^" 1 + 2.52 x 10- 2 T 3 - 2 ) 

(Jj-L 

- 4.2 x 10- 9 s- 1 ^ 5 ^ 4 exp(^^) [2.5 x (IQ^ 1 + lO' 2 ^ 2 )] , (6) 

where we adopted C p = 1.4 x 10 8 erg g -1 K -1 , T 3 is the temperature T in units of 
1000 K, Vfas is the zonal velocity scale in units of km s _1 and B\ is the magnetic field 
strength in G. 

The scaling in Eq. (JBJ) reveals that instability is favored at lower atmospheric densi- 
ties (hence lower pressures), for larger values of the magnetic field strength and the zonal 
wind velocity scale, and at higher temperatures (T > 10 3 K), with a strong exponen- 
tial T-dependence. For conditions representative of the mbar pressure level on the day- 
side of a hot Jupiter such as HD209458b, T ~ 1500 K, p ~ 10 -8 g cm -3 , one obtains 
dH ohm /dT ~ 2.9 x 10~ 3 s _1 for B = 10 G and V$ = 3.3 km s -1 . By comparison, the radia- 
tive relaxation time r ra d ~3x 10 3 s at the mbar level in such an atmosphere (e.g., Showman 
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et al. 2008), which indicates thermally unstable conditions since dH ^ m /dT > l/r ra d- The in- 
stability would disappear, however, for lower values of the magnetic field strength (B < 1 G), 
the wind velocity scale (V^ < 1 km s _1 ) or the atmospheric temperature (T < 1000 K). 

This derivation is useful to illustrate the potential for thermal instability and to un- 
derstand parameter dependencies but it is only valid when magnetic drag is inefficient and 
when thermal ionization is dominated by potassium. In order to establish more reliable 
domains of instability, we now relax these assumptions and turn to a numerical parameter 
space exploration. 



3. Thermo- Resistive Instability Domain 

Method - Relaxing the weak drag limit implies the need for a magnetic drag law that 
specifies the magnitude of drag as a function of atmospheric resistivity, r\ and magnetic field 
strength, B. In the interest of simplicity, we adopt here a parameterized drag law inspired 
from the physical arguments developed in Menou (2012). We evaluate the dragged zonal 
velocity scale, V£, from the drag-free value according to 

yd = Y± ( 7 \ 

where Vd ra g = -Rp/^drag is a characteristic drag velocity evaluated from the drag time 
Tdrag (previously defined in Eq. [2]) and the planetary radius, R p . We use the dragged 
velocity, Vjf, in all our ohmic dissipation calculations. This formulation satisfies the weak 
drag expectation — > at low temperatures (high resistivities) and the strong drag 
exponential dropoff, V? oc r\ (for n e — 1), expected at high temperatures (e.g Menou 2012). 
The additional dimensionless parameters n s and n e , with default values of unity, are used to 
control the threshold and the steepness of the transition from the weak to the strong drag 
regimes. With this drag law, the half-speed transition for a P = 60 mbar pressure level (near 
the thermal photosphere) occurs at a temperature T ~ 2250 K for B = 3 G, T ~ 1900 K for 
B = 10 G and T ~ 1650 K for B = 30 G, with a half-width AT ~ 200-300 K. Adopting a 
weaker drag law, with n s = 0.3 and n e = 1, pushes these transitions to higher temperatures 
(T ~ 2500, 2000 and 1750 K for B = 3, 10 and 30 G). While this parameterized magnetic 
drag law is broadly consistent with the results discussed in Menou (2012), it also clearly is 
one of our most uncertain model ingredients. 

We use here a more detailed, numerical solution to the Saha equation for the thermal 
ionization fraction x e (described in Menou 2012), rather than the analytic, potassium-only 
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formulation discussed in the previous section. This allows us to explore more consistently the 
low resistivity and high temperature (T > 1700 K) regime which is relevant to the hottest 
exoplanets. Finally, we adopt a zonal velocity profile, with default scale — 7 km/s (before 
drag is applied), that decreases exponentially with depth according to exp(— Pi), where Pi 
is the pressure in units of 1 bar. This velocity scale and its rapid decline at high pressures 
are broadly consistent with the results of Rauscher & Menou (2012) and other published 
circulation models on the properties of superrotating equatorial winds on hot Jupiters. 

We identify instability domains by numerically evaluating the ohmic derivative term 
dHohm/dT and comparing it to the inverse of the radiative relaxation time, r ra d- Regions 
satisfying the inequality in Eq. (T5|) are considered unstable. We estimate the radiative 
times on the basis of the work of Iro et al. (2005) and Showman et al. (2008), which 
yield consistent results for HD209458b. We note that the tabulated, planet-specific values 
provided by Showman et al. (2008) approximately follow the simple scaling r ra d oc T~ 3 
over a wide range of pressures, as may be expected from the formal definition of r ra d (e.g., 
Goody & Yung 1989). In the interest of generality, we normalize our radiative times to the 
values found by Iro et al. (2005), using fits provided in Heng et al. (2011), and at a each 
pressure level, we scale the radiative time in proportion to (T iro /T) 3 , where T iro (P) is the 
corresponding temperature value in the profile of Iro et al. (2005). While this procedure is 
not expected to yield r ra d values more accurate than a factor of a few, it is sufficient given 
other sources of uncertainty in our modeling methodology. 

Results - Figured] presents dayside temperature-pressure profiles for a prototypical hot 
Jupiter with radius R p = 10 10 cm and surface gravity g = 900 cm s -2 . Fifty logarithmically- 
spaced profiles are shown for values of the planetary effective temperature from T e g = 1000 K 
to 2500 K. The profiles are calculated with the one- dimensional radiative solution presented 
by Guillot (2010), for a zero incidence angle and a dilution factor / = 0.5 representing a 
dayside average. An internal heat flux corresponding to T int = 150 K is adopted. As shown 
by Guillot (2010), even when adopting constant thermal and visible absorption coefficients 
as is done here (ft t h = 10~ 2 cm 2 g _1 and k v = 4x 10~ 3 cm 2 g" 1 ), these temperature-pressure 
profiles are reasonably good solutions for hot Jupiter atmospheres over a large range of 
insolation fluxes. For reference, HD209458b has T e g ~ 1420 K which corresponds to the 
twentieth profile from the left. 

Thermo-resistive instability domains are also shown in color in Fig. [1] for a default 
magnetic drag law (n s = n e = 1 in Eq. [3]). The green and red instability domains correspond 
to a zonal velocity scale = 7 km s _1 and magnetic field strengths of B = 10 and 40 G, 
respectively. The blue instability domain corresponds to V$ = 10 km s" 1 and B = 3 G. 
For V<j, = 7 km s -1 and B = 3 G no instability is found. As expected, higher velocities 
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promote instability. Larger magnetic field strengths also promote instability and tend to 
shift the instability domain to lower atmospheric temperatures. Details of the adopted T-P 
profiles (e.g., opacity coefficients, dilution factor /) have only minor effects on the instability 
domains in the sense that instability is primarily set by physical conditions which correspond 
to specific regions of the T-P plane. 

The boundaries of the instability domains, which can be extended in pressure but are 
typically limited to a narrow range of temperatures, can be understood as follows. To 
the left of each instability domain, at low temperatures, resistivity varies too weakly with 
temperature for instability to occur. At the bottom of each instability domain, besides the 
rapid exponential drop of velocity with pressure, we find that our assumption of a radiatively- 
dominated atmosphere breaks down. We thus arbitrarily truncated the instability domain 
at the point where r a d v = Rp/V^ = r rad , with the implicit assumption that strong advective 
cooling would likely neutralize the instability. More stringent criterions (r ra d <C T a d v ) would 
further reduce the instability domains to pressure < 0.3-0.1 bar in all our figures. Finally, 
to the right and at the top of each instability domain, the emergence of the strong magnetic 
drag regime rapidly leads to V£ <C and the effective disappearance of the instability. 

Figure |2] shows that a weaker drag law (with n s = 0.3, n e = 1 in Eq. [3]), for the same 
zonal velocity scale = 7 km/s and magnetic field strength B = 10 G, can result in a 
significantly more extended instability domain (compare to the green domain in Fig. [TJ), by 
pushing the transition to strong drag to higher temperatures. Using n e > 1 in Eq. (E]) to 
steepen the transition to strong drag has only minor effects on the instability domains. 

Interestingly, the upper atmosphere of hot Jupiters and other strongly-forced exoplanets 
are subject to strong photo-ionization rates. An examination of the arguments in section 
§0 shows that the addition of a temperature-independent contribution to the ionization 
fraction does not impact the value of dH ^ m /dT or the instability criterion. However, the 
lower resistivity due to the additional ionization results in stronger magnetic drag, which is 
a limiting factor for the instability. Figure [3] shows the same instability domain as in Fig. |2j 
when a constant ionization fraction = 10~ 8 is added everywhere in the atmosphere, to 
emulate the effects of photoionization. This level of additional ionization stabilizes the lowest- 
density, upper atmospheric regions, where magnetic drag is the strongest. We find that 
the entire atmospheric domain can be stabilized with the addition of a constant ionization 
fraction x+ ~ 10~ 7 . 

Fortney et al. (2003) modeled sodium photoionization in the atmosphere of HD209458b 
and found ionization fractions in excess of thermal values at pressure levels < 30-100 mbar, 
corresponding to x+ > 10~ 8 (for a solar abundance of sodium). This, together with the 
results shown in Fig. [3J suggests that plausible levels of photo-ionization may indeed con- 
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trol the occurrence of the thermo-resistive instability in hot exoplanet atmospheres. At 
particularly high photo- ionization levels (and possibly super-solar atmospheric abundances), 
photo-ionization could entirely suppress the instability Provided one associates the thermo- 
resistive instability with the observationally-inferred phenomenon of dayside temperature 
inversions, these results would appear consistent, at least qualitatively, with the observa- 
tional trend established by Knutson et al. (2010): hot Jupiters subject to the strongest 
photo-ionization fluxes are the least likely to show such inversions. 

4. Discussion and Conclusion 

Our argument for the existence of a thermo-resistive instability in the atmospheres of 
hot, strongly-forced exoplanets relies on a number of assumptions which are not all well 
justified. Most importantly, the induction formalism used here and in other studies of mag- 
netic effects in hot Jupiter atmospheres and interiors assumes axisymmetry and steadiness. 
We already emphasized in Perna et al. (2010a) the clearly non-axisymmetric properties of 
strongly-forced exoplanet atmospheres, a situation that would only be aggravated by the 
onset of a localized thermal instability (e.g., on the dayside). 

Our study of an MHD-based instability with a steady-state induction equation is another 
important shortcoming. Dimensional analysis of the time-dependent axisymmetric induction 
equation (see Eq [1] in Menou 2012) suggests that induced currents would lag changes in 
temperature and resistivity by a resistive diffusion time, Tdm ~ L 2 /i]. Adopting a pressure 
scaleheight for the characteristic lengthscale L, we estimate that r diff < r rad for the green 
and red instability domains in Fig. [TJ which suggests that current-adjustment delays will not 
strongly affect the instability development. At high temperatures and low densities, however, 
around the blue instability domain and at higher temperatures in Fig. [H resistivities become 
small and r^g > r ra d, which may invalidate the implicit assumption of instantaneous current 
adjustment made in our steady-state evaluation of the ohmic dissipation term in Eq. fl2]). 
A careful consideration of this issue likely requires a full MHD treatment and we shall 
simply note here that time- dependent current adjustments may impact the thermo-resistive 
instability of the hottest exoplanet atmospheres beyond the simple treatment adopted here. 

Despite such limitations, it is tempting to associate the thermo-resistive instability 
mechanism identified here with the dayside thermal inversions inferred for a number of 
hot Jupiters. As shown in Fig. (TJ the instability domain is restricted to a specific range of 
pressures and temperatures. For magnetic field strengths B ~ 3-30 G, instability domains 
can match well the upper atmospheric dayside conditions of moderately hot exoplanets such 
as HD209458b, with low enough nightside temperatures T < 1100 K to suppress the insta- 
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bility. On cooler hot Jupiters such as HD 189733b, dayside conditions would barely achieve 
instability (see T-P profiles in Showman et al. 2008 and Rauscher & Menou 2012). By 
contrast, very hot exoplanets with T c g > 1500-2000 K may be too hot on their daysides 
for instability, and even if their nightsides were to meet the temperature requirements, the 
instability may be suppressed because of the very effective wind drag exerted on the dayside 
or possibly strong deviations from radiative equilibrium on the nightside. In principle, some 
hot exoplanets may be preferentially unstable near their limb, which will be an interesting 
issue to explore with improved instability models. 

The outcome of the instability is difficult to anticipate beyond qualitative expectations. 
Unstable regions should reach temperatures high enough for saturation of the instability 
in the strong drag regime, which can amount to temperature excesses of several hundreds 
Kelvin according to the width of domains shown in Figs. dH3J Perhaps more importantly, the 
radiative response of the vertically-coupled atmospheric layers, together with the horizontal 
coupling caused by locally modified thermal and drag conditions, will result in a very non- 
linear response that is best studied with global models. While short instability growth times 
may prove numerically challenging, global circulation models with adequate MHD treatments 
offer a promising avenue for progress. They would help clarify the global energetics of the 
instability, which is ultimately powered by the same differential thermal forcing as the global 
circulation itself. 

The author thanks E. Rauscher for comments on the manuscript. This work was sup- 
ported in part by NASA grant PATM NNX11AD65G. 
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Fig. 1. — Thermo- resistive instability domains superposed on temperature-pressure profiles 
representative of the dayside of a typical hot Jupiter, for the default magnetic drag law 
(n s = n e = 1 in Eq. [3]). The green and red instability domains correspond to a zonal 
velocity scale = 7 km s _1 and magnetic field strengths B = 10 and 40 G, respectively. 
The blue instability domain corresponds to V<f> = 10 km s _1 and B = 3 G. 
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Fig. 2. — Same as Figure 1 with a weaker magnetic drag law (n s = 0.3, n e = 1 in Eq. [3]), 
= 7 km s -1 and B = 10 G. This drag law delays the transition into the strong drag 
regime to higher temperatures, which extends the domain of instability (compare to the 
green domain in Figure 1). 
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Fig. 3. — Same as Figure 2 with a constant ionization fraction = 10 added throughout 
the atmosphere, to mimick the effects of photoionization. The additional, constant ionization 
contribution stabilizes the upper atmospheric regions by inducing a stronger drag regime. 
The entire atmospheric domain can be stabilized with the addition of a constant ionization 
fraction x+ ~ 10~ 7 . 



